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ABSTRACT 

It is generally accepted that gamma-ray bursts (GRBs) are initiated by a 
relativistic pair fireball, converting its internal energy into kinetic energy of a 
relativistically moving plasmoid and subsequently into radiation. Here, we inves- 
tigate the early stages of this evolution, after the pair fireball has become optically 
thin to 77 pair production. We show that for a short period of time, ~ 0.1 - 
a few seconds after the initial explosion, the pair plasmoid evolution might be 
dominated by coUisional processes prior to the formation of a coUisionless shock. 
We simulate these processes during the early pair plasmoid evolution and cal- 
culate the expected radiative signatures. We show that the radiation from the 
collision-dominated pair plasmoid phase results in a short (~ a few ms) flash of 
thermal soft X-ray emission, followed by a transition phase of <,1 s during which 
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the fireball turns Thomson thin, but its radiation remains dominated by thermal 
Comptonization, peaking at around E'pk ~ 100 MeV - a few GeV. While the 
very early thermal emission could be associated with the quasi-thermal radiation 
signatures found in the very early phases of several bright BATSE GRBs, the 
predicted subsequent flash of high-energy emission should be easily detectable 
with the GLAST satellite. 

Subject headings: gamma rays: bursts — gamma-rays: theory 



1. Introduction 

With the estabhshment of the cosmological distance scale of 7-ray bursts (GRBs) and 
the considerable success of the synchrotron-shock model (Meszaros & Rees 1993; Katz 
1994; Tavani 1996) to explain the broadband (radio through X-ray) continuum afterglow 
signatures of GRBs (Wijers, Rees, & Meszaros 1997; Vietri 1997; Waxman 1997; Galama 
et al. 1998; Meszaros, Rees, & Wijers 1998), it is now generally accepted that GRBs are 
initiated by a relativistic pair fireball, transferring its internal energy into kinetic energy of 
a relativistic blast wave (Cavallo & Rees 1978; Shemi & Piran 1990; Meszaros & Rees 
1992; Rees & Meszaros 1992). Non-thermal particle acceleration at the forward shock 
and subsequent radiative cooling is believed to lead to the formation of a broken power- 
law distribution of ultrarelativistic electrons, producing the observed broadband afterglow 
radiation primarily through synchrotron emission (e.g.. Sari, Piran, & Narayan (1998); 
Bottcher & Dermer (2000)), though Compton upscattering of synchrotron photons may 
also play a significant role (Dermer, Bottcher, & Chiang 2000; Sari & Esin 2001). 

In the relativistic fireball/blastwave model for GRBs, a large amount of energy is released 
in a short period of time. This initial configuration is highly opaque to 77 absorption and 
pair production and leads to the formation of a relativistic pair fireball. As the fireball 
expands, it cools adiabatically, converting its internal energy into bulk kinetic energy of the 
outflow, and its 77 opacity decreases until it becomes optically thin to 77 pair production. 
At this time, the pairs are essentially cold in a reference frame co-moving with a small section 
of the expanding pair plasma / radiation shell. It has been recognized very carly-on (e.g., 
Shemi & Piran (1990); Meszaros & Rees (1992)) that even a small intrinsic or external 
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contamination of the pair fireball with baryons would lead to re-conversion of most of the 
fireball energy into kinetic energy of the outflow, and that any radiation escaping the pair 
fireball during this phase of energy conversion would have quasi-thermal signatures. 

However, more detailed measurements of time-averaged photon spectra emerging from 
GRBs, in particular by the Burst and Transient Source Experiment (BATSE) on board 
the Compton Gamma-Ray Observatory (CGRO) have shown that they can generally be 
fitted successfully with a model consisting of two power-laws, smoothly connected by an 
exponential turn-over (Band et al. 1993). Such a power-law high-energy spectrum was 
interpreted as the radiative signature of non-thermal particle acceleration at relativistic 
shocks formed in the process of deceleration of the relativistic blast wave produced during 
the initial fireball stage (e.g., Paczyhski & Rhoads (1993); Meszaros & Rees (1993)). 
While it is now generally agreed that the smoothly decaying afterglow emission observed 
from many GRBs is produced by relativistic electrons accelerated at the external forward 
shock, in the stage when the blast wave is efficiently being decelerated by the sweeping-up of 
external material, many researchers now believe that the emission during the prompt GRB 
phase may be produced in internal shocks during collisions of subsequent relativistic blast 
waves produced in a series of energy release events from the central engine (e.g., Fenimore, 
Madras, & Nayakshin (1996); Sari & Piran (1997)), although an external-shock scenario for 
the prompt GRB phase cannot be ruled out at this time (e.g., Dermer & Mitman (1999)). 

Both internal and external shock scenarios for GRBs are generally starting out with 
the assumption that the relativistic blast wave has produced a relativistic shock with fully 
developed hydromagnetic turbulence in order to allow for efficient non-thermal particle ac- 
celeration at the shock front. This neglects the fact that it takes a finite amount of time 
to build up the necessary turbulence in the two-stream multi-fluid system consisting of the 
relativistically moving pair plasma and the background hydrogen plasma. While a detailed 
treatment of the process of the development of hydromagnetic turbulence through a rela- 
tivistic two-stream instability is rather cumbersome (see, e.g., Pohl & Schlickeiser (2000); 
Schlickeiser & Dermer (2000); Schlickeiser et al. (2001)), the time scale for the develop- 
ment of such turbulence can be estimated by the Alfven crossing time through the shell of 
relativistically moving pair plasma. Using hydrodynamics simulations, Meszaros, Laruna, & 
Rees (1993) have demonstrated that a relativistic pair flreball becomes rapidly compressed 
into a shell of co-moving thickness A' with R/A' ~ Fq during the early evolution of the 
flreball, and that the ratio R/A' remains approximately constant after this compression 
phase, as long as the blast wave is coasting (i.e. F const.). A relativistic pair flreball of 
this thickness, carrying a total energy per unit solid angle of Eq = 10^^ -Eq 52 erg sr"-*^ = 
meC^n'p^^^R^ = 8.2 x lO^^n'i^Rl^ erg sr~^ is initially highly optically thick to Thomson 
scattering and thus radiatively inefficient since radiation will remain trapped within the 
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shell. Here, n'^g^^,. — 10^^ cm~^ is the co-moving pair density (throughout the paper, 
primed quantities refer to the reference frame co- moving with the pair plasmoid), and 
i?i5 = 10^^-Ri5 cm, where ~ 10^^ cm is a typical value for the Thomson thinning radius (see 
§6). We assume that the magnetic field given by the equipartition parameter = 0.1 6^,-1 
between magnetic field and pair plasma energy density (i.e., u'^ = esUp^^j.), and that the 
pair plasma is cold in the co-moving frame. Then, for a bulk Lorentz factor F = 10^ Fa, the 
Alfven crossing time scale in the observer's frame is 



^ 3 ^B,-l 

Thus, Eq. (1) indicates that it may take ~ 0.1 - 1 s, i. e. a non-negligible fraction of 
the prompt GRB phase, to build up the hydromagnetic turbulence for non-thermal particle 
acceleration at either internal or external shocks associated with these pair blastwaves. 

The time scale (in the observer's frame) for coUisional processes, tcoib be estimated 

as 



tcou ~ (Fcnp^i.cTT) ^ ~ 5 X 10 * (n'14 F3) ^ s, (2) 

indicating that during the initial phase — prior to the generation of strong hydromagnetic 
turbulence — the radiative and energy- exchange processes in the blast wave are dominated 
by coUisional processes. This fact is particularly interesting in view of the evidence in a 
significant fraction of GRB that the early time-resolved BATSE spectra show too hard a 
low-energy spectrum to be produced by non-thermal synchrotron emission (Crider et al. 
1997; Preece et al. 1998). Recently, Preece (2001) has found evidence in several BATSE 
GRBs that their very early spectra are described better by thermal rather than nonthermal 
emission. 

Setting both time scales (1) and (2) equal, one can estimate the condition under which 
non-thermal particle acceleration may begin to dominate the particle energization behind 
the forward shock. Wwe find that coUisional processes will dominate as long as 
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Although alternative emission models for the prompt GRB radiation have been proposed 
by several authors (e.g., Brainerd (1994); Liang (1997)), the dynamics and radiative signa- 
tures of the transition phase between the time when the fireball becomes optically thin and 
the time of efficient non-thermal particle acceleration by hydromagnetic turbulence has not 
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been studied in detail. In this paper, we are simulating the evolution of and time-dependent 
emission from a relativistic pair plasmoid before the generation of hydromagnetic turbulence 
leading to non-thermal particle acceleration. In §2, we describe and motivate the basic as- 
sumptions of our model. The computation of the evolution of the thermal pair plasma and 
the fate of the swept-up background plasma are outlined in §3. In §4 we describe the calcu- 
lation of the time-dependent radiation spectra emitted by the resulting particle distributions 
in the blast wave. The evaluation of the global dynamics of the blast wave is given in §5. 
Results of our numerical study are presented in §6. We summarize in §7. 



We assume that the pair plasmoid forming the blastwave is a circular shell segment 
with opening solid angle fl. Initially, it consists of a pure pair plasma with a pair density 
'^pair ~ ^^^^ ^^14 cm~^ in the co-moving frame of the pair plasmoid. At the time when the 
pair fireball becomes optically thin to 77 pair production, the pairs are essentially cold in the 
co-moving frame. The blastwave interacts with an external medium (normal electron-proton 
plasma) of density nisM = 1 ^t-ism.o cm~^ in the stationary frame of the surrounding medium. 
At any given time, the plasmoid is located at a radius R — 10^^ i?i5 cm from the center of the 
explosion and has a thickness A' = R/T in the co-moving frame. The bulk Lorentz factor of 
the plasmoid is denoted by F = 10^ T-^. The background magnetic field in the blastwave is 
B = 1 Bq G and parametrized by the equipartition parameter 63 = £5,-1 through 



where ©pair = kT^air/ {fTieC?) is the normalized (co-moving) pair temperature in the plasmoid, 
and Kn is the modified Bessel function of 2nd kind (a.k.a. McDonald function) of order n. 
For simplicity, we assume eb to be constant throughout our simulation. In fact, a gradual 
build-up of the magnetic field from lower initial values of may extend the period of 
applicability of our approach. 

We are following the evolution of both the pair plasma as it is being energized by 
Coulomb and inelastic (i.e. bremsstrahlung) collisions with swept-up electrons and protons 
from the surrounding medium, and cools via radiative and adiabatic cooling. Simultaneously, 
we solve for the evolution of the population of suprathermal protons and electrons which are 
swept up from the environment with an initial Lorentz factor F in the co-moving frame of 
the plasmoid, and transfer their energy to the background pair plasma and into radiation. 



2. Model setup 
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3. Evolution of pairs, electrons and protons in the blastwave 

Electrons and protons which are swept up by the pair blastwave have an initial Lorentz 
factor r in the co-moving frame, and will transfer part of their kinetic energy to the back- 
ground pairs in the plasmoid via elastic scattering (M0ller/Bhabha and Coulomb scattering, 
respectively), thus heating the thermal pair plasma. The remainder of the kinetic elec- 
tron/proton energy will be radiated away through bremsstrahlung with the background 
plasma and through synchrotron radiation. The thermal background plasma will be heated 
by collisions with the swept-up electrons and protons and cool via bremsstrahlung and syn- 
chrotron emission, adiabatic losses and pair annihilation. The very short time scale for 
elastic scattering in the plasmoid (see Eq. [2]) indicates that the incoming particles are 
rapidly isotropized in the co-moving frame of the plasmoid. Thus, we assume local isotropy 
of the particle distributions. 

3.1. Suprathermal electron/proton cooling rates 

In the following, we quote simple approximation formulae for the cooling rates for 
suprathermal electrons and protons in a cold background pair plasma, which will be used in 
our simulations to follow the evolution of the swept-up non-thermal particle distributions. 
For highly relativistic test particles, the elastic scattering cross sections for electron-electron 
(M0ller) and electron-positron (Bhabha) scattering are approximately equal, and we may 
use a simple power-law fit to the M0ller scattering energy exchange rate of a test particle 
with a cold background pair plasma given by Nayakshin & Meha (1998). We find 

7M0,e ~ 2.6 X 103^4 7,-°-°2s-\ (5) 
The e~e~and e~e+bremsstrahlung energy loss rate can be approximated as 

7br,e~8xlO-'n;4 7,^-^^s-^ (6) 
Suprathermal electrons will also emit synchrotron radiation, thus losing energy at a rate 

Finally, we take into account synchrotron-self- Compton emission and estimate the associated 
energy loss as 
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where u'^^ ^ L'^^ ^/ [Vt c) is the energy density in synchrotron photons. These estimates 
indicate that except for the highest-energy electrons (with 7e^lO^), the energy loss will 
be dominated by elastic scattering, leading to rapid thermalization of the suprathermal 
electrons. 

For Coulomb losses of suprathermal protons in the cold pair plasma, we use Eq. (4.22) 
of Mannheim & Schlickeiser (1994) with j3 ^ j3(. (i.e. ultrarelativistic protons in a non- 
relativistic or mildly relativistic background plasma): 

7coui,p ~ 3.3 X 10-2 ni4S-i (9) 

The bremsstrahlung energy loss rate of suprathermal protons is evaluated by integrating 
the respective photon spectrum calculated by citetjonesTl over the outgoing photon energies. 
Since to our knowledge, the result has not yet been published elsewhere, we derive the 
complete expression in Appendix A. In our simulations, we use a simple power-law fit to the 
resulting bremsstrahlung energy loss rate: 

7b,,p« 1.1 X 10-5 n;4 7j=^s-i (10) 

Finally, we take into account synchrotron emission of ultrarelativistic protons, leading to an 
energy loss of 

As in the case of suprathermal electrons, we expect that the swept-up protons will lose energy 
primarily via elastic scattering with the thermal background pair plasma. 

3.2. Suprathermal electron and proton spectra 

As mentioned in the previous section, we generally expect that the energy loss of both 
protons and electrons might be dominated by Coulomb interactions for 7p,max = Tr^^lO^- 
However, for the highest-energy protons and electrons, synchrotron losses may also play an 
important role. The evolution of the suprathermal electron and proton populations may be 
calculated by solving the continuity equation for the cooling particles: 
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- {iN[j,t']) + NU^,t'), (12) 

where Ns„{'~f,t') = /^rcF njsM^ -R^ ^(7 — T) is the sweep-up rate of ISM particles, and we 
have neglected escape. Since the cooling time scale for ultrarelativistic particles is of the 
order of or shorter than the dynamical time scale, we may approximate the evolution by a 
sequence of quasi-equilibrium solutions to Eq. (12), given by 

/ N 0.02 

(^) for 1 < 7e < 7e''^ 

Nehe) =N'^{ (13) 

for j^' < 7e < r 



, 2 

7^' 



for the electrons and 



for 1 < 7p < 7p'" 

Npi^p) =K{ , ^ -2 (14) 

(^) for 7^^"^ < 7p < r 

for the protons. The break energies are defined through the condition 7eiast.scat. = 7sy and 
given by 

7^'" = 8.5 X 10^ {n[4)^ (15) 

and 

7?"^ = 3.95 xlO«(n;jV2 5-1. (16) 

If 7gyp > r, the high-energy branches of Eqs. (13) and (14) are not realized, and the 
respective particle spectra are single truncated power-laws for 1 < 7e/p < F. 

The background pair plasma is assumed to maintain a thermal distribution at a tem- 
perature determined by the balance between heating by the swept-up particles and radiative 
cooling via bremsstrahlung and synchrotron emission. 



3.3. Temperature of thermal pair plasma 

As mentioned above, most of the kinetic energy of the protons is transferred to electrons 
through Coulomb collisions with the thermal background pair plasma. Thus, the heating 
rate of the thermal plasma is given by 



-9- 



Lheat ^^R^r^ rupc^ pcnisM ~ 4.50 x 10^^ fl Rl^ Vl nisM,o erg 8"^ (17) 
This heating should be balanced by cooling through bremsstrahlung, 



L 



br,th 



3.42 X 104« Q {n'^^f Rl, ./e~ erg s'^ for Opair < 1 

1.36 X 10^9 Q (^/ j2 ^/^ Q^^.^ (ln[2epair] + 0.673) erg s'^ for Op^ir > 1 



(18) 

— where A'12 = A'/(10^^ cm) — , synchrotron emission, 

L,y,th ^ 3.175 X 10^1 n n[, R\, A[, Bl 6 erg s"^ (19) 

(cf. Bottcher et al. (1999)), and thermal Comptonization of the soft photon field (at 
E <^ rrieC^). For typical parameters, the pair plasmoid is expected to be moderately optically 
thick, tt^I, and have mildly relativistic temperatures, ©<,1. In this regime, we may use 
the analytic solution of Hua & Titarchuk (1995) for the saturated-Comptonization case 
(their Eq. [9]). We assume that the dominant soft input radiation field is the electron 
synchrotron component, which we approximate, for simplicity, by a 5 function spectrum 
Ls(e) = i^sy,e '^(c~f^^c )• The cooling rate due to Thermal Comptonization is then calculated 
by integrating the Comptonized spectrum (in the comoving frame) over photon energy. 

Adiabatic cooling is described by a virtual luminosity Ladi = 3£'thc/i?, where E'th 
is the thermal energy content of the thermal pair plasma. At nonrelativistic pair plasma 
temperatures, pair annihilation becomes important. We use the expressions of Svensson 
(1982) for the pair annihilation cooling rate and the pair annihilation rate to calculate the 
cooling and depletion of the number of pairs due to this process. 



4. Radiation spectra 

The radiation emitted by the thermal pair plasma can be calculated using standard 
expressions (see, e.g., Bottcher et al. (1999), and Hua & Titarchuk (1995) for the Thermal 
Comptonization component). The intrinsic luminosities of the bremsstrahlung and synchro- 
ton components are given by Eqs. (18) and (19), respectively. 

The (optically thin) bremsstrahlung spectrum of the suprathermal electrons in the co- 
moving frame reproduces the broken power-law shape of the electron spectrum: 
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for 6' < et,^, 

^br,e(^') = ^^r,e <! , , , -1 (20) 

{£) for<,<.'<r 

where e^J! g = 7^"" and e' = hu' /{rrieC^). The suprathermal proton bremsstrahlung spectrum 
may be approximated as 

r 1 for e' < i 

^;r,p(^) = Ll,, (21) 

U2eO-^-'-^ fori<e'<r 

(Jones 1971). 

The nonthermal synchrotron spectra can be approximated as 



T' (f'\ = 

-^sy,e/pl"^ I -^sy,e/p 



1/3 

for e' < e^' 



sy,e/p 



sy>e/p _ 

1/2 

for e^'' ; < e' < e°^^^ 



where 



with = 4.414 X lO^'^ G. 



(22) 



C = I- (7e'0' (23) 



Cy!? = £ (24) 

br 3 rUe B / br\2 /qcN 

'-p=2;;^5:^^^) ^^^^ 



Using a simple delta-function approximation for the Compton-scattering cross section 
in the Thomson limit, the evaluation of the electron-SSC spectrum is straightforward, but a 
bit lengthy. The full expression is given in Appendix B. 

The normalization factors of the radiation spectra are given by setting the energy- 
integrated luminosities. 
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T _ rO ^br 



1 + ln 



r 

^br,e 



(27) 



^br,p 



(i + 4[i-(2r)-o-25]) 



,br 



-br 



'sy,e/p 



rbr 

'sy,e/p 



^sy,e/p -^sy,e/p ^ g ''sy.e/p ^ V "=sy,e/p 

-^SSC.e = -^SSC ^ 

where N is given in Appendix B (Eq. [B4]), equal to the fraction of swept-up kinetic particle 
energy, which is transferred into the respective radiative cooling channel, 



(28) 

(29) 
(30) 



2 2 2 lO O'br c 

Lbr,e = r i? rUeC (3 CUiSM -. —. "^r —. 

7M0,e + 7br,e + 7sy,e + 7SSC,e 


(31) 


T o T-i2 d2 2 a 7br,p 

Lbr.p = ytV R rupc pcriisM . . . 

7Coul,p + 7br,p + 7sy,p 


(32) 


L,y,e ^nr^R' rriec' 13 c msu -. ^ ■ ^1^'! — — 

7M0,e + 7br,e + 7sy,e + 7SSC,e 


(33) 


7Coul,p ~r 7br,p + 7sy,p 


(34) 


T m-'2 d2 2 a 7SSC,e 

7M0,e + 7br,c + 7sy,c + 7SSC,e 


(35) 



All luminosities have been calculated in the co-moving frame and are transformed to apparent 
isotropic luminosities through 



1 — pvn '''77 {^') 

i^app(6) = h ^ L'{e'), (36) 

wherein = (r[l — f3rcos6ohs])~^ is the Doppler beaming factor, e = De', and = min(l, ct'/A')x 
min(l, l/[r^r2]) is a correction factor accounting for the following two effects: First, due to 
light-travel-time constraints, only emission from a fraction of ~ ct' /A' of the plasmoid will 
be visible to the observer. Second, if the opening angle of the cone in which plasmoid is 
moving, is larger than the beaming angle l/F, then only a fraction l/(r^Q) of the surface 
will contribute to the observed emission. In this analysis, we neglect the fact that in reality 
the observed emission is a superposition of emission from different layers of the plasmoid 
which have different evolutionary ages. This effect may alter the detailed spectral shape 
and light curves slightly (see, e.g., Panaitescu & Meszaros (1998) for a detailed analysis of 



- 12 - 



this effect in the case of GRB afterglows), but will not change the quahtative results of our 
analysis. 

The effect of pair production by 77 absorption is treated by injecting a number of 
thermal pairs equal to the number of absorbed photon pairs within any given time step 
into the plasmoid. Since we have shown that thermalization is the dominant energy-loss 
process for suprathermal particles, this approach yields an accurate description of the 77 
pair production process in the plasmoid. 



5. Plasmoid deceleration 

Due to momentum conservation, the bulk Lorentz factor F of the pair plasmoid will 
decrease as it sweeps up external matter, according to 

dr = (37) 

(Dermer & Humi 2001) where dm = i?^ (mp + nie) nisu P c dt, rfi^adi is the loss of internal 
energy due to adiabatic cooling, and M is the total, relativistic mass in the fireball (i. e. 
rest mass + internal kinetic energy of pairs + mass of swept-up external material). The 
rest-mass increment dm is related to the increment in M by 

dM^rdm-^^^^^±^ (38) 
where (iE'rad is the net energy produced in radiation throughout the plasmoid. 



6. Numerical results 

For a general parameter study, we have performed a series of pair plasmoid simulations 
with various values of bulk Lorentz factors F, co-moving pair densities n'pg^^^, initial radii Rq, 
and external matter densities riisM- Throughout our simulations, we have fixed the magnetic 
field equipartition parameter to es = 0.1, and the observer is assumed to be located on the 
symmetry axis, i.e. ^obs — 0. We keep the ratio R/A' — Fq constant at its initial value. 

As briefiy mentioned in the introduction, our simulations generally start out in a highly 
Thomson thick regime, in which the output is basically a thermal blackbody at the blue- 
shifted pair temperature. For this reason, the results are virtually independent of the initial 
radius Rq, and for the results presented in the following, we have chosen Rq — 10^^ cm. 
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In this initial phase, the plasmoid is inefficient in terms of high-energy emission. At the 
radius where the Thomson depth is ~ 1, the plasmoid becomes radiatively efficient, and a 
flash of high-energy 7-ray emission is produced. We denote the time at which the plasmoid 
becomes marginally Thomson-thin as tx- This transition time is plotted as a function of the 
essential model parameters -^^,52 and Tq in Fig. 1 which shows that it gradually increases 
with increasing injected pair energy and decreasing bulk Lorentz factor. We find typical 
values of 0.01 s <tT<0.1 s. 

For each simulation, we then evaluate, among other quantities, the apparent quasi- 
isotropic uLi, peak luminosity and the peak photon energy £'pk at the time of maximum 
received flux. The dependence of those quantities on the parameters £'^,52 and Fq is illus- 
trated in Figs. 2 and 3. 

In these calculations, we have fixed nisM = 100 cm^^. The graphs show that we find 
typical peak i/Li, luminosities ~ 10"'^ - lO^'' ergs s^\ increasing with increasing Ea,52 and 
r. uFi, peak energies E'pk are typicaUy in the range ~ 100 MeV - several GeV, and increase 
with r, while they are only weakly dependent on the injected energy Eq^^2- 

For a typical simulation in this range of parameters. Fig. 4 shows the composite spec- 
trum with the individual radiation components at the time of maximum observed fiux. The 
X-ray and 7-ray spectrum is clearly dominated by thermal Comptonization in this phase, 
time corresponding to Fig. 4, the radial Thomson depth of the plasmoid is tt ~ 0.24; the 
pair temperature is ©pair ~ 0.44. The spiky shape of the electron and proton synchrotron 
spectra is due to the 5 function approximation to the respective emissivities. 

Fig. 5 illustrates the spectral evolution during the coUision dominated phase of the 
pair plasmoid evolution. It shows that the hard X-ray and 7-ray spectral output remains 
dominated by thermal Comptonization throughout most of this phase. During the first few 
ms, there is a strong soft X-ray component due to the thermal emission of the cold pair 
plasmoid, which is rapidly shifting to higher energies as the pairs are being heated. At very 
high energies, E>^1 GeV, there is a weak tail from suprathermal proton bremsstrahlung. 
According to Eq. (1), hydromagnetic turbulence might develop within ~ 1 s in the case 
discussed here. 

Light curves expected from the collision-dominated pair plasmoid are shown in Fig. 7. 
The figure illustrates that the maximum spectral power output is expected around ~ 0.1 s 
after the formation of the pair plasmoid. Generally, the time of maximum fiux decreases 
with increasing photon energy, indicating an overall hard-to-soft spectral evolution. 
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7. Summciry and conclusions 

We have investigated the early stages of the evolution of a relativistic pair fireball, im- 
mediately after it has become optically thin to 77 pair production. We have shown that 
for a short period of time, ~ 0.1 - a few seconds after the initial explosion, the pair blast 
wave evolution might be dominated by coUisional processes prior to the formation of a col- 
lisionless shock which would subsequently lead to non-thermal particle acceleration at the 
shock front. We have simulated the relevant energy exchange and radiation processes, in- 
cluding suprathermal electron and proton thermalization and bremsstrahlung, suprathcrmal 
electron and proton synchrotron emission, thermal pair bremsstrahlung, thermal pair anni- 
hilation, thermal Comptonization, and 77 pair production and absorption, during the early 
pair plasmoid evolution, and calculated the expected radiative signatures. 

We have investigated the dependence of the radiative signatures on pair plasmoid pa- 
rameters which are very hard to predict from first principles in the framework of current 
GRB progenitor models. In the range of total injected energies of Eq ~ 10^^ ergs sr^^ and 
initial bulk Lorentz factors To ~ 500 - 1000, the radiation from the collision-dominated pair 
plasmoid phase results in a short period (~ a few ms) of thermal soft X-ray emission from 
the initially Thomson thick, cold pair plasmoid. As the pair plasma is heated due to the 
sweeping-up of external material and is becoming Thomson thin due to expansion, the ob- 
servable emission turns into a quasi-thermal spectrum, peaking at around E'pk ~ 100 MeV - 
a few GeV, dominated by thermal Comptonization by the mildly relativistic pair plasma in 
the plasmoid during the first ~ 0.01 - 1 sec after the onset of the GRB. The apparent peak 
pLi, luminosities are expected in the range of ~ 10^^ - 10^*^ ergs/s, sustained over typically 
a few tens of milliseconds. 

The expected very early, thermal signatures may already have been observed in time- 
resolved, early BATSE spectra of a few GRBs. The quasi-thermal spectra expected during 
the optical-thinning transition phase, ~ 0.01 -Is after the onset of the GRB, are peaking 
in the ~ 100 MeV regime. Thus, they may have remained undetectable for BATSE, but 
should be easily detectable by the GLAST mission, scheduled for launch in 2005. 
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A. Bremsstrahlung losses of suprathermal protons 

In order to calculate the energy loss rate of suprathermel electrons due to bremsstrahlung 
in inelastic collisions with a cold pair plasma, we start out with the differential cross section 
^ as given in Eq. (4) of Jones (1971) (for consistency we use e for the normalized energy 
of the outgoing photon, which is a in Jones' notation) and evaluate 

^max 

I ^ rrif, [ , da , . ^ , 

7br,p = -n^ PpC — / dee— (Al) 
TTir) J de 



Up 





where 



e^ax = 1.123 7p/3|e--. (A2) 
Evaluating the integral in (Al), we find 

where o; ~ 1/137 is the fine structure constant, and 

GH„) = - (\nAf+ -ln2H In^H tt^ H ln2 + - , (A4) 

^'^^ 3^ ^ V3 18/ 27 9 40 \A J ' ^ ^ 

where A 0.687^ ^ 1 for relativistic protons. The last term in Eq. (A4) denotes all terms 
of order l/A <C 1. 



B. The electron-SSC spectrum 



We evaluate the electron-SSC spectrum using a S function approximation for the Comp- 
ton scattering cross section in the Thomson regime (Reynolds 1982), with a sharp cut-off 
at incident photon energies ejg > 1. In the case 7^'' = T this yields: 



In the case 7g'' < F, we find: 



3 Bcr 



-Drr 
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(B2) 



The integrated luminosity in this spectrum is 



-^ssc 



(B3) 



where 
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Fig. 1. — Thomson thinning time of a relativistic pair plasmoid injected at Rq = 10^^ cm, 
as a function of total injected energy per unit sohd angle, E^, and initial bulk Lorentz factor 
Fq. Other parameters: riisu = 100, and R/ A' = F. The green dots indicate simulated values; 
the surface has been constructed using a spline interpolation. The surface colors encode the 
values of (vertical axis in the lower panel) with red corresponding to ~ 0.5 - 1 s and 
blue corresponding to ~ 0.5 - 1 ms. 
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Fig. 2. — Peak vLy apparent isotropic luminosity at the time of maximum received flux, as 
a function of £"^,52 and Fq. Parameters are the same as in Fig. 1. The green dots indicate 
simulated values. The surface colors encode the values of z/Lj^ peak with red corresponding to 
T-'L„^peak ~ 10^^ ergs s~^ and blue corresponding to vLy^^eak ~ 10^^ ergs s~^. 
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Fig. 3. — Peak photon energy at the time of maximum received flux, as a function of -En,52 
and Tq. Parameters are the same as in Fig. 1. The green dots indicate simulated values. 
The surface colors encode the values of -Epeak with red corresponding to -Epeak ~ 10 GeV and 
blue corresponding to -Epeak ~ 100 MeV. 
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Fig. 4. — Composite photon spectrum from the blast wave at the time of maximum received 
flux (tree = 64 ms). Parameters for this simulation are: Fq = 10^, Up^ir ~ x 10^^ cm~^, 
Rq — 10^^ cm, implying -Eq,52 = 1; nisM = 100 cm~^, = 0.1, 6'obs = 0. 77 absorption has 
been taken into account to calculate the total emission, while the individual contributions are 
plotted without correction for 77 absorption. Thus, the importance of 77 absorption and pair 
production is illustrated by the difference between the individual, unabsorbed components 
and the total, emerging spectrum at E>,1 GeV. 
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Fig. 5. — Time evolution of the observed photon spectra for the same set of parameters as 
in Fig. 4. 
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Fig. 6. — Time evolution of several relevant quantities describing the evolution of the plas- 
moid, for the same set of parameters as in Fig. 4. 
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Fig. 7. — Light curves at several different observed plioton energies for tlie same set of 
parameters as in Fig. 4. 



